Infectious Disease Modelling

Brazil experienced its worst dengue epidemic in 2024, with an increase in severe outcomes among elderly individuals, highlighting the need for vaccination strategies informed by population immunity. We reconstructed state- and age-specific immune profiles by fitting an age-structured catalytic model to dengue surveillance data from 2000-2025. This inferred immune profile was incorporated into an SIR-type transmission model to project disease burden and to evaluate the public health impact of vaccine allocation strategies. We observed substantial spatial heterogeneity in accumulated immunity, which was closely associated with projected burden and optimal allocation. Strategies prioritizing states with lower historical incidence and larger susceptible populations yielded greater benefits. Under limited vaccine supply, targeting children resulted in the greatest reductions in DALYs, whereas targeting elderly populations achieved smaller DALY reductions but larger reductions in mortality. We further showed that accelerated vaccine supply markedly increased impact, while allocation choices remained critical when availability was constrained.

Brazil has one of the highest dengue burdens globally, with transmission influenced by climate events. Yet, the effects of compound climate extremes on dengue risk and prediction have not been quantified to date. Using a Bayesian spatiotemporal model of monthly dengue cases across 435 health regions from 2010 to 2025, we estimated relative risk associated with interacting long-lag and short-lag drought indicators and temperature anomalies across nine Koppen-Geiger climatic zones, adjusting for absolute temperature and El Nino-Southern Oscillation. Long-lag dry, mid-lag warm and short-lag wet conditions were associated with the highest relative risk in five climatic zones, with temperature anomalies modifying hydrometeorological dynamics more strongly in temperate than tropical areas. Climate covariates improved the prediction of case counts and outbreaks compared to baseline models at lead times of 1, 3, 6 and 12 months, and for the upcoming season 5 months ahead. Improvements were greatest in the Northeast, Southeast, and Midwest. Recent dengue incidence provided the greatest improvement at the 1-month lead, whereas combining incidence with climate information performed best at 3- and 6-month leads. We extended our framework to weekly state-level predictions of the season ahead to align with national public health planning. We proposed a multi-staged dengue forecasting approach to enhance early warning and response across Brazil.

Introduction Despite the significant burden of Shigellosis, routine surveillance systems in resource-limited settings lack predictive capabilities to guide proactive public health responses. Consequently, we developed and evaluated the accuracy, uncertainty and computational cost of multi-horizon forecasting models for Shigella incidence among children <5 years in rural western Kenya using 16-year longitudinal incidence time-series. Methods We estimated monthly Shigella incidence using data from three pediatric enteric studies (2008-2024), supplemented with Kenya Health Information System during inter-study gaps. Incidence estimates incorporated sampling, healthcare-seeking, stool-testing, and pathogen-specific weights, with bootstrap-derived 95% confidence intervals. Predictors included Google Trends health-related indicators, climate, temporal trends, and seasonality. Forecasts at 1-, 3-, and 6-month horizons were evaluated using rolling-origin validation. Following horizon-specific predictor selection using Boruta, 14 state-space, machine-learning, and hybrid models were compared using point-forecast and probabilistic metrics, bootstrap confidence intervals, and Model Confidence Sets. Adaptive conformal prediction recalibrated prediction intervals, while Shapley Additive exPlanations assessed predictor contributions. Results Mean monthly Shigella incidence was 2.51 per 100 child-years (standard deviation=1.85). Shigella dynamics in Kenya can be forecast beyond seasonal persistence alone for up to six months reducing absolute errors by ~ 33%. However, optimal model architectures and predictor sets shift systematically across forecasting horizons. Temporal trend and recent incidence were the dominant predictors, while predictor importance shifted with forecast horizon toward syndromic, pharmaceutical, seasonal, and lagged environmental signals. Extreme Gradient Boosting performed best at 1 month (Root mean squared error [RMSE]=1.86), remained competitive at 3 months (RMSE=1.89), whereas exponential smoothing (ETS) performed best at 6 months (RMSE=2.05). Forecast accuracy declined with increasing horizon. ETS and TBATS provided the best-calibrated prediction intervals, whereas machine-learning models generally under-covered. Model Confidence Sets demonstrated horizon-dependent model distinguishability. Conclusion Horizon-specific model selection balancing predictive accuracy, uncertainty calibration, and computational feasibility could strengthen operational early warning systems for shigellosis.

We develop an age- and immunity-structured COVID-19 model with setting-specific contacts, calibrated to Ontario incidence and hospital admissions before Omicron. We simulate counterfactual vaccination strategies and contact restrictions in schools, workplaces, and other community settings, evaluating projected severe-case demand against an illustrative healthcare-capacity benchmark. Vaccination combined with workplace and community contact restrictions keeps simulated demand within capacity. In contrast, relaxing restrictions in either setting produces large surges and repeated capacity threats, even when school restrictions remain. School-only changes produce smaller effects on severe disease. These findings highlight the value of vaccination and setting-specific contact restrictions that prioritize workplace and community settings for pathogens with adult-concentrated severe outcomes.

Wolbachia-based interventions were less effective in achieving introgression in Rio de Janeiro than other settings. Reduced introgression occurred in locations where dengue was historically high, raising concerns of inequity. We examined factors affecting Wolbachia introgression in Rio through trajectory modeling, spatiotemporal Bayesian models, and Ae. aegypti abundance estimation. During the intervention (2017-2020), Wolbachia introgression was 26% (median; IQR,14-44%), with four distinct trajectories. Despite receiving more releases, densely populated areas with low socioeconomic status had lower introgression. In spatiotemporally explicit multivariable regressions, higher Ae. aegypti abundance, temperature, population density and lower socioeconomic status were associated with reduced introgression. Improved characterization of Ae. aegypti abundance and other structural/environmental predictors of introgression could be used to calibrate release doses to ensure equitable benefits from Wolbachia interventions.

Japanese encephalitis virus (JEV) is a leading cause of viral encephalitis in humans in Asia. Prior to 2021, only sporadic outbreaks had been detected in Australia's tropical northern extremes. However, during 2021-23, JEV infection was identified in 45 people and in more than 85 piggeries, together with many detections in feral pigs, mosquitoes, and other livestock across Australia, including in the temperate far south. Here we estimated the spatial extent of virus transmission during the outbreak (the realised niche) and the area environmentally suitable for transmission (the potential niche). We collated a national dataset of JEV surveillance and detection in humans, mosquitoes, domestic pigs, feral pigs, and sentinel chickens from 2021-23. We combined these data with (1) previously published predictions of habitat suitability for vertebrate host and vector species, and (2) estimated densities of humans, domestic pigs, and feral pigs, in a multipart assessment of spatial risk. Our model projected broad potential environmental suitability over much of northern and eastern Australia during the outbreak period, but floodplain wetland areas in the north and the inland east of Australia had both the highest estimated environmental suitability for transmission, and the strongest evidence for virus circulation.

Importance: Coccidioidomycosis, also known as Valley fever, is a serious fungal disease caused by inhaling Coccidioides spores aerosolized from the soil. Despite rising public health concern, the distribution and incidence of coccidioidomycosis in the US remain poorly characterized as not all states are required to report cases, and case under-ascertainment is presumed to be high, particularly outside of the endemic region. Objective: To estimate monthly, county-level coccidioidomycosis incidence across the US while accounting for widespread under-ascertainment. Design, Setting, and Participants: We integrated multisource data, including human cases diagnosed in electronic health records, serological testing of dogs, and ecological and sociodemographic predictors, in a machine learning model estimating monthly, county-level coccidioidomycosis incidence. The model was trained on surveillance data from Arizona and California from 2013 to 2023, then applied to predict incidence in states across the contiguous US over the same time period. Main Outcomes and Measures: Incidence of coccidioidomycosis at the monthly county level. Results: The model predicted coccidioidomycosis incidence from 2013 to 2023 for 3,108 counties from all states in the contiguous US. The model identified 60 counties in non-reporting states where predicted incidence was [&ge;]2 cases per 100,000 population (considered moderate to high incidence). In states with mandatory reporting outside of Arizona and California, the model predicted county-level incidence rates that were 28% higher (95% bootstrap interval, 0.1%-69.8%) than those captured in national surveillance, equating to a total of 608 potentially missed cases per year (95% bootstrap interval, 340-964). Conclusions and Relevance: In this cross-sectional study, the model quantified case under-ascertainment for coccidioidomycosis relative to Arizona and California, and highlighted the potential for alternative data sources to fill in gaps in mandatory reporting. Estimating incidence while accounting for under-ascertainment and non-reporting may help efforts to improve public and physician awareness and to detect emergence in new regions.

Tick-borne disease risk reflects the combined influence of vector abundance, pathogen prevalence, human exposure, and seasonal activity, yet these components are rarely integrated at the county level. We developed species-specific ecological risk scores for three medically important tick species (Acari: Ixodidae): Amblyomma americanum (L.), Dermacentor variabilis (Say), and Ixodes scapularis Say, across all 102 Illinois counties using data from 2019 to 2022. County-level predicted tick abundance and pathogen prevalence were scaled and combined using equal weights to generate primary vector-level risk scores. Robustness to the equal-weight assumption was evaluated by systematically varying the relative weights assigned to abundance and pathogen prevalence, while a separate sensitivity analysis incorporated reported tick-bite incidence where available. Associations with human disease were evaluated using negative binomial regression of tick-borne disease case counts on the corresponding risk scores, with county population included as an offset. Monthly phenology indicators were constructed by integrating human disease cases, reported tick bites, and sampled tick abundance and were evaluated in relation to temperature and precipitation. Seasonal risk was estimated by combining county-level vector risk with species-specific monthly phenology, generating monthly risk surfaces that were deployed through an interactive R Shiny application. Higher risk scores were significantly associated with greater disease case rates for each tick species. County rankings were highly stable under moderate changes in abundance and prevalence weighting. Seasonal activity differed among species, with A. americanum and D. variabilis showing concentrated warm-season risk and I. scapularis exhibiting a broader seasonal profile. These county-level seasonal risk maps identify where and when relative tick-borne disease risk is elevated in Illinois, supporting targeted prevention, surveillance, and clinical awareness.

To deal with the COVID-19 pandemic, Denmark implemented an aggressive program aimed at controlling the pandemic which included both testing and contact tracing programs. The large amount of data collected in Denmark on the pandemic and its integrated nature made it possible to investigate the effectiveness of the testing and contact tracing effort. The results of that investigation and a discussion of the effectiveness and limitations of the different control efforts are reported herein, as are lessons learned concerning possible improvements for the future. The analysis is limited to the period from the start of the COVID-19 pandemic in Denmark (late Feb. 2020) into May 2021, at which time the impact of vaccination intervention was already playing a role with ~25% of the Danish population at least partially vaccinated. From the available data sets we estimate the national Infection Detection Efficiency (IDE) starting around July 1, 2020, the time at which the national infection detection program started to be implemented in a more robust form. The results rely on daily COVID-19 hospital admissions, PCR test results, and anti-body prevalence values for the Danish population. We find that for the long-term average, the combined Danish national testing and contact tracing program identified on average about 55 to 60% of the infected individuals during the period from Sept. 1, 2020, into early May 2021. However, higher infection detection efficiencies were achieved in localized areas for short time intervals. Our investigations further indicate that the Danish testing and contact tracing program, which was among the most ambitious in the world, could have had a better balance between voluntary probing on the one hand and focused contact tracing and targeted testing on the other hand. As typically implemented, it proved to be largely ineffective in controlling the pandemic. However, in at least one case when it was implemented in an enhanced and focused mode, it was able to bring a large local outbreak that dominated the country's cases under control. We have compiled a list of lessons learned regarding the use of contact tracing procedures to control potential future pandemics before vaccinations become available.

Background Influenza vaccines have seen many changes over recent years in terms of both policy recommendations and product types. Immune responses and vaccine effectiveness vary by vaccine strain composition and product type. This study aimed to explore this variation on a regional level through examining spatial and temporal autocorrelations of vaccine uptake of different product types within Michigan from 2017 to 2021. Methods Influenza vaccine administration data from the 2017/18 to 2021/22 seasons were obtained from the Michigan Care Improvement Registry. Population size estimates for 2017 and 2019 were obtained from the Michigan Department of Health and Human Services. Age-standardized and age-stratified vaccination rates were computed. Disparate county-level vaccine uptake clusters were constructed using spatial autocorrelation indicators and spatiotemporal clusters were identified using scan statistics. Linear mixed effect models were used to study the impact of social vulnerability and urbanicity level on seasonal influenza vaccination rate. Results Vaccination incidence rates for adjuvanted, recombinant, and nasal sprays vaccines and for any vaccines increased significantly through the study period. County-level spatiotemporal clusters of type-specific uptake rates were identified. Most notably, we identified an emerging cluster with high uptake of adjuvanted vaccine radiating from Ottawa County, located in suburban Grand Rapids, and the disappearance of high uptake cluster for high-dose vaccine uptake in the area. Vaccination rates by any product type were found to vary seasonally and showed weak negative association with more rural regions. The overall vaccination rates appear to be lower among counties with higher social vulnerability in certain urbanicity levels, but social vulnerability and urbanicity were not significantly associated with overall vaccination rate during the study period. Conclusions Vaccine type-specific uptake was spatially and temporally diverse in Michigan between 2017 and 2022. Mid and Northern Michigan had consistently relatively lower vaccination rate compared to other regions in Michigan.

Background. Nigeria's reported Lassa fever case fatality ratio (CFR) rose from 18.9% at week 23 of 2025 to 25.0% in 2026, prompting concerns over failing clinical care. Because surveillance artefacts can also alter CFRs, candidate mechanisms that published data can resolve were identified and evaluated. Methods. This ecological study extracted weekly state-level counts from 474 Nigeria Centre for Disease Control and Prevention situation reports (2017-2026). Per a pre-specified plan, years were compared at phase-matched points, year-on-year CFR changes were decomposed into within-state and between-state components, and tipping-point ascertainment shifts calculated. Weekly incident and cumulative published totals were compared post hoc. Findings. Of 11 candidate mechanisms, 2 could be evaluated fully and 5 in part. On weekly incident counts, 2026 did not differ from 2025 at the phase-matched midpoint (+0.3 percentage points, 95% CI -5.0 to 7.6); within-state (-1.1, 95% CI -4.7 to 3.0) and between-state (+0.5, -2.0 to 2.9) components spanned zero. Conversely, the published cumulative CFR was 5.4 to 6.2 percentage points higher because 70 deaths were added retrospectively by week 36 outside weekly reports. Reconciling cumulative figures would require 279 missed survivors (6% of suspected cases not confirmed). Interpretation. The 2026 CFR rise came from deaths added to cumulative totals outside the weekly counts, and published data cannot attribute it to care. Endemic zoonotic CFRs cannot serve as clinical outcome measures unless testing volumes, deaths added late, and post-mortem sampling are published alongside case counts. Funding. US NSF-NIH-NIFA EEID; UK BBSRC; Birgitta Sintring Fellowship (Uppsala University).

Routine infant vaccination has helped reduce cases and deaths from diphtheria, tetanus, and pertussis globally. However, these diseases continue to pose serious health risks as protective immunity wanes over time. Three booster doses of diphtheria-tetanus-pertussis-containing vaccine (DTPCV) administered from 12-months to 15-years of age increase duration of protection beyond childhood. Growing demands on government health spending underscore the importance of understanding the health and economic impact of vaccines prior to introduction. Using an integrated multi-disease epidemiological-economic model, we estimated the health impact (cases averted, deaths averted, and changes to population immunity) and budget impact of introducing DTPCV booster doses in Uganda. Compared to no DTPCV boosters, introducing all three booster doses is projected to result in 2.2 times the proportion of the population protected against symptomatic diphtheria and 2.8 times against tetanus infection among Ugandans under 15 years of age by 2039. Compared to introduction of the early childhood booster (EC) at 18 months of age, introducing all three DTPCV boosters simultaneously is projected to avert 3.2 times the number of tetanus cases over fifteen years at 2.4 times the incremental cost. The projected government vaccination cost per planned dose was USD 0.58 for simultaneous introduction of all three boosters compared with USD 0.68 for the EC booster alone. The average annual incremental cost of all three DTPCV booster doses is USD 2.79 million. Under modelled assumptions, simultaneous introduction of all three DTPCV booster doses in Uganda provides the greatest health gains, with a lower government vaccination cost per planned dose than introduction of the EC booster alone, but may still pose affordability challenges for Ugandas government.

Vegetation plays a critical role in shaping fine-scale tick density by structuring local microclimates, habitat resources, and host interactions that largely determine tick survival and activity, with implications for tick-borne disease transmission. Tick burden in the United States continues to increase, with annual reported tick-borne disease cases more than doubling from 2004 through 2019. These trends call for a broader scope of tick-mitigation strategies that go beyond traditional prevention measures that focus on personal protection; and there is an opportunity to inform habitat-level interventions. Here, we quantify fine-scale vegetation structure across 48 subplots at tick surveillance sites in Western Massachusetts to assess local habitat features associated with Ixodes scapularis density and potential for human-tick contact risk. In assessing understory and overstory vegetation attributes, this study identified life stage-specific differences in predicted tick density, with higher adult densities in areas with deeper leaf litter and presence of invasive plant species and nymphs associated with high presence of ferns. Both adult and nymphal Ix. scapularis were predicted at higher densities in relatively low basal area forests (18.9 to 26.9 m2/hectare). This study design provides a simple, low-cost framework for organizations looking to incorporate tick-informed forest measurement ecology into their public health tick surveillance protocols with minimal training.

Introduction: Raised blood pressure is an important cardiovascular risk factor in Cambodia; however, national averages may obscure substantial local geographic heterogeneity. We assessed the spatial distribution of raised blood pressure and its demographic, metabolic, behavioural, built-environment, and environmental correlates using the 2023 Cambodia STEPwise approach to noncommunicable disease (NCD) risk factor surveillance (STEPS) survey. Methods: We analysed 289 georeferenced primary sampling units (PSUs) from the 2023 Cambodia STEPS Survey. PSU-level case counts and numbers of adults with valid blood pressure data were modelled using a frequentist binomial generalised additive model with a two-dimensional Gaussian-process smooth. An explanatory model included selected PSU-level and contextual covariates; a separate mapping model used spatially continuous predictors to generate predictions on a 5-km grid. Design-based national and residence-specific prevalence estimates were calculated separately using the survey design and Step-2 weights. Results: Among 4,225 adults with valid blood pressure data, unweighted prevalence of raised blood pressure was 28.6%, compared with 18.8% after Step-2 weighting; weighted prevalence was 17.5% in rural areas and 21.7% in urban/suburban areas. Mean observed unweighted PSU-level prevalence was 29.0%. In the adjusted spatial explanatory model, the proportion aged 60-69 years ({beta}=1.969, SE=0.327; p<0.001) and overweight/obesity prevalence ({beta}=0.996, SE=0.265; p<0.001) were positive ecological correlates. The spatial model had lower Akaike information criterion and residual deviance than its non-spatial counterpart, although Bayesian information criterion was higher. Prediction uncertainty was greater in some sparsely sampled areas. Conclusions: Raised blood pressure showed significant geographic heterogeneity across Cambodia. Older age composition and overweight/obesity prevalence were positive PSU-level correlates, not individual-level causal effects. Predicted prevalence should be interpreted alongside uncertainty and should complement, rather than replace, design-based STEPS prevalence estimates for geographically informed NCD planning.

BACKGROUND Timely cause-of-death (COD) information guides health policy, yet many deaths in low- and middle-income countries are not medically certified. In Africa, fewer than 10% of deaths are. Verbal autopsy (VA) yields a structured symptom questionnaire and a free-text narrative. Mortality surveillance increasingly records the narrative alongside the questionnaire, but the algorithms used to produce VA-based statistics were designed for structured responses only, so narratives inform physician COD assignment but do not contribute to algorithmic COD assignment. We evaluated whether supervised models could use both modalities to improve COD assignment. METHODS We analyzed 4,999 VA records for deaths at ages 12 years or older from a national sample of registered deaths in South Africa (RSA), using physician-assigned CODs as the reference standard. We trained unimodal and multimodal pipelines on the labelled records: pretrained language models (PLMs) fine-tuned on narratives alone, machine learning (ML) models on structured responses alone, and the two combined at the data, feature, and decision levels. We also prompted a locally deployed generative large language model (LLM) on the same unimodal and multimodal inputs, without parameter updates. All were compared with InSilicoVA, the operational comparator representing current automated practice. We evaluated individual-level assignment, population-level agreement, calibration, and generalization across provinces in RSA, reporting these metrics rather than ranking pipelines on any one of them. RESULTS Every supervised pipeline improved on InSilicoVA by 8.5-15.9 percentage points in top-1 accuracy and 3.5-7.1 points in cause-specific mortality fraction (CSMF) accuracy. The narrative and structured pipelines performed similarly in aggregate but differed by cause: 63.9% of deaths received identical predictions, 11.6% were correct only from the narrative pipeline, and 13.6% only from the questionnaire pipeline, with the narrative pipeline more accurate for noncommunicable causes and the questionnaire pipeline for communicable causes. Against the unimodal approaches, the decision-fusion pipelines added 4.2-5.4 percentage points in top-1 accuracy, concentrated in the most common causes, whereas feature fusion with cross-modal attention improved the less common causes, adding 6.6 points (95% CI, 1.7-11.5) in balanced accuracy and 5.5 points (95% CI, 0.6-10.3) in macro F1. Feature fusion reached the highest CSMF accuracy (95.7%). Paired differences from the better unimodal representative ranged from -0.2 to +3.4 points, with all intervals including zero. The LLM exceeded InSilicoVA on individual-level metrics and reached the highest balanced accuracy of any pipeline (60.6%), while the supervised pipelines retained higher CSMF accuracy. Supervised pipelines had lower expected calibration error than InSilicoVA. Supervised gains held across all nine held-out provinces, including every top-1 comparison and nearly every macro F1 and CSMF comparison. CONCLUSIONS PLMs using narratives alone substantially improved on current best practice and matched ML pipelines using structured responses, while multiple multimodal strategies produced further gains in individual COD assignment. The zero-shot LLM showed particular strength on rare causes, whereas supervised pipelines better recovered the population cause distribution. These findings support incorporating narratives as a core input to automated VA systems and evaluating representative fusion strategies for the intended output. External validation across settings and age groups will support broader implementation.

Background: TB vaccine development is hampered by the natural history of M. tuberculosis: most infected people rarely progress to disease, so trials must be very large and long to accrue sufficient disease endpoints. Current trials enrol IGRA-positive participants (IGRA: Interferon-Gamma Release Assays), but IGRAs remain positive for many years and cannot distinguish recent from remote infection. Tests for recent infection are now in development, and we argue they could make trials substantially more efficient. Here, we investigate how screening participants prior to enrolment with a test for recent M.tb infection could impact the total number of enrollees required for an adequately-powered TB vaccine trial. Methods: Using a statistical model of vaccine trial recruitment, we estimated the number needed to contact, test, and enrol to achieve an 'effective' trial (defined as 50 disease endpoints), under IGRA-gated enrolment (life-long duration of positivity) vs. enrolment with a test for recent infection. We considered five recruitment strategies: non-test-based enrolment for ages 15-49, ages 12-34, and household contacts of recent cases; and test-based enrolment using IGRA, and a test for recent infection with two-year duration of positivity. The model explored trial recruitment in South Africa, Kenya, and India. Results: Compared with an extended-positivity test, a test for recent infection reduced the median number needed to enrol by 4.2- to 6-fold under favourable conditions. It also outperformed non-test-based strategies, including enrolling household contacts. This comes at the cost of substantially more testing, but crude estimates suggest recency-based enrolment would still be cheaper overall provided the new test does not cost more than roughly 5-10 times an IGRA. Conclusions: Diagnostics for recent M.tb infection have the potential to improve vaccine trial efficiency, but more research is now needed to understand test characteristics and to evaluate their use in this setting. Our Viewpoint sets out the test attributes that would need to be established, such as time to positivity, duration of positivity, sensitivity, and specificity, before such tests could be used in vaccine trials.

Background: Both over- and under-antibiotic use pose major global health challenges, yet there is no consensus on how to quantify deviations from optimal prescribing or set evidence-based targets. We developed a benchmarking tool for inpatient antibiotic use and evaluated its application in 10 countries. Methods: We developed decision-tree models to estimate expected empiric antibiotic use for 11 indications, based on the WHO AWaRe book for two age groups: 28 days to 12 years and >12 years. We parameterised the models using country-specific antimicrobial resistance estimates in the 10 countries and cross-sectional Global Point Prevalence Survey (Global-PPS) case-mix data from 2021 for 10 countries and from 2015-2022 in a subset of five countries. We then compared estimated expected antibiotic use with observed empiric use (from Global-PPS) in defined daily doses for adults and days of therapy for children, across indications and AWaRe categories. Findings: In 2021, among patients aged >12 years, in all countries but one, observed empiric Access antibiotic use, as a proportion of total antibiotic use, was lower than expected and empiric Watch antibiotic use was higher than expected. Similar patterns were observed in all countries but one in patients aged 28 days to 12 years. The magnitude of discrepancies between expected and observed Access and Watch antibiotic use varied across indications and countries. The median observed-to-expected ratios of usage across countries ranged from 0.3-0.8 for Access antibiotics, and 0.7-5.3 for Watch antibiotics in patients aged >12 years in 2021. Temporal trends from repeated surveys from 2015 to 2022 showed little change in observed-to-expected ratios for Access and Watch antibiotics. Interpretation: This benchmarking tool uses Global-PPS data to quantify differences between observed antibiotic use and WHO AWaRe recommendations. Our tool can be adaptable to local guidelines, provides policy-relevant estimates of expected antibiotic use to support monitoring, procurement planning, and targeted stewardship intervention.

Background: Severe malaria remains a major cause of hospitalization in malaria-endemic sub-Saharan Africa, and anemia is among its most frequent severity criteria. Data on the impact of anemia on hospital length of stay (LOS) among adults with severe malaria remain limited in Central Africa. We assessed the impact of anemia on LOS among adults hospitalized for severe malaria in Brazzaville, Republic of Congo. Methods: We conducted a cross-sectional, descriptive, and analytic study of medical records from adult patients ([&ge;]18 years) diagnosed with severe malaria (WHO criteria) and admitted to four hospitals in Brazzaville between January 2020 and December 2024. We dichotomized LOS as short (<5 days) or prolonged ([&ge;]5 days). Multivariable logistic regression, using a backward stepwise procedure, identified factors independently associated with prolonged LOS; isolated hemoglobinuria served as the reference severity category. Results: Among 523 patients hospitalized for severe malaria, 296 met the inclusion criteria. Median age was 25.5 years (IQR 20-40); 51.7% were female, and 82.1% had a low socioeconomic status. Median hemoglobin was 10.55 g/dL, and median LOS was 4 days (IQR 3-5). In multivariable analysis, prolonged LOS was independently associated with delay to consultation (aOR 1.25/day, 95% CI 1.09-1.44), isolated anemia (aOR 7.07, 1.23-40.55), convulsions (aOR 5.86, 1.81-18.97), prostration (aOR 4.42, 1.21-16.2), and impaired consciousness (aOR 4.29, 1.17-15.7). The strongest associations involved combined severity criteria: anemia with prostration (aOR 8.58), anemia with jaundice or hemoglobinuria (aOR 16.21), and impaired consciousness with anemia (aOR 33.31). Conclusions: Anemia is a major, potentially modifiable determinant of prolonged hospitalization in adults with severe malaria, with a markedly amplified effect when combined with neurologic impairment. Early anemia screening, improved access to safe blood transfusion, and reduced delays to care could shorten hospital stays and ease the burden on health systems in this setting.

Background: Filoviruses, specifically Ebola virus (EBOV), Bundibugyo virus (BDBV), and Marburg virus (MARV), have been the etiologic agents of multiple outbreaks of severe disease in the Democratic Republic of the Congo (DRC). Two EBOV-specific vaccine regimens have been approved for use in outbreak scenarios, including the Ad26.ZEBOV/MVA-BN-Filo (Zabdeno/Mvabea) two-dose vaccine, where the MVA-BN-Filo booster dose contains EBOV, MARV, and Sudan (SUDV) virus glycoproteins (GPs) and the nucleoprotein from Tai Forest virus. While long-term immunogenicity of up to two years has been found in clinical trial settings, limited real-world data exists. Methods: Using a multiplexed immunoassay, samples (n=698 at baseline) collected in Mbandaka, Equateur province, DRC, between September 2023 and October 2024 were assessed for seroreactivity to six filovirus antigens (1 BDBV, 2 EBOV, 2 MARV, and 1 SUDV) at three distinct periods: baseline/pre-vaccination, 6-8 weeks post-first dose (Zabdeno), at booster administration (Mvabea), and one-year post-boost dose. Results: Substantial filovirus seroreactivity was found at baseline, with 48.7% (n=340/698) of individuals being seroreactive to at least one virus antigen. High baseline seroreactivity was found against BDBV GP (31.1%, n=217/698), with lower baseline seroreactivity to EBOV GP (15.5%, n=106/698), comparatively. Overall EBOV GP seroreactivity increased to 66.8% (n=376/563) at one year-post vaccine boost and 57.6% (n=250/434) of participants were seroreactive. Despite the booster dose inclusion of SUDV GP and MARV GP, there were no substantial increases in reactivity against either antigen at subsequent timepoints. Overall filovirus seroreactivity increased at all timepoints for the entire cohort, although filovirus reactivity decreased over the study period for those who were seroreactive at baseline. Findings: These data highlight detection of significantly elevated seoreactivity to EBOV GP following first-dose administration of Ad26.ZEBOV, which remained relatively stable one-year post-booster dose administration. Additionally, this study identifies high levels of baseline seroreactivity to BDBV GP despite no previous documented outbreaks in this region of DRC. These data suggest that vaccination with the Ad26.ZEBOV-MVA-BN-Filo vaccine induces stable seroreactivity against EBOV GP, but not to other filovirus antigens.

Background: Treponema pallidum quantitative polymerase chain reaction (qPCR) positivity and cycle threshold (Ct) by rapid plasma reagin (RPR), disease stage, and other characteristics are poorly characterized in routine-care specimens, leaving it unclear when qPCR of such specimens is informative for diagnosis. Methods: This Japanese multicenter cross-sectional study explored qPCR positivity and Ct by these characteristics. Nonblood specimens collected during routine care underwent qPCR at a central laboratory for this study. Bayesian logistic and normal regressions modeled positivity and mean Ct of TP0574 and polA with stage-specific continuous automated RPR terms, patient characteristics, collection site, and facility random intercepts. Estimates were standardized over the analyzed specimens. Population-group predictions were instead averaged within each group. Results: Among 154 nonblood specimens from 145 patients (71.7% non-MSM men and women) at 14 facilities, 109 (70.8%) were qPCR positive. At RPR 1, adjusted positivity was 72.2% (95% credible interval, 60.9%--82.9%) in primary and 43.8% (16.7%--69.9%) in secondary syphilis. Positivity showed a modest downward trend with RPR in primary but increased in secondary syphilis, reflecting a stage-by-RPR interaction (secondary-to-primary ratio of odds ratios, 1.78; 1.25--2.65). Predicted positivity showed an upward trend with RPR values among MSM; adjusted positivity was relatively high for genital/urethral/urinary and anal/perianal/rectal specimens. Adjusted mean Ct increased with RPR in primary and followed an inverted-U pattern in secondary syphilis. Conclusions: In primary syphilis, nonblood qPCR positivity was high even at low RPR. When serologic findings remain inconclusive in this stage, qPCR is highly recommended for securing an early diagnosis.